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Abstract 

Ab initio computations within the density functional theory are 
reported for the spin cross-over complex, [Fe(btz) 2 (NCS) 2 ], where 3d 6 
Fe 7/ is characterized by High Spin (HS t 2ff , e^) and Low Spin (LS 
t 2s , e°) states. Results of infra red and Raman spectra for the iso- 
lated molecule are complemented for the crystalline solid with a full 
account of the electronic band structure properties: the density of 
states assessing the crystal field effects and the chemical bonding as- 
signing a specific role to the Fe-N interaction within the coordination 
sphere of Fe^ . 

1 Introduction 

Few inorganic transition metal ion complexes have the possibility of exhibit- 
ing two electronic states of d electrons, the High Spin (HS) and Low Spin 
(LS). Switching between these two states is subjected to small energy mag- 
nitudes around ~ \lbt and the transitions can be achieved with external 
constraints such as temperature, pressure as well as by applying light. Such 
spin cross-over (SCO) behaviour has been discussed in detail in several re- 
views [H EJ [3j E] . Theoretical studies at the molecular level [SI E] are equally 
available but to our knowledge there are no solid state investigations of the 
electronic structure for the whole crystal system. Such studies are likely to 
provide further information on the changes of the electronic density of states 
as well as on the chemical bonding in the close neighbourhood of central 
Fe ion. The purpose of this work is to provide spectroscopic and electronic 
structure information on SCO system [Fe(btz) 2 (NCS) 2 ] (btz stands for 2.2'- 
Bi-4.5-dihydrothiazine), both for the isolated molecule and for the extended 
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solid state. The btz ligand is similar in its structure to 2,2'-bipyridine-N,N' 
aromatic ligand. However the presence of sulfur replacing carbon provides 
a non-planar bidentate ligand as it is shown by the projection in fig. [TJ 
Lastly NCS~ mono-dentate ligands are in a cis-position like in most of SCO 
complexes, such as [Fe(phen) 2 (NCS)2] [7j. 

The crystallographic and SCO characteristics were investigated by many 
authors [HI El LTQ] . Just like the prototype of SCO systems, [Fe(phen) 2 (NCS)2], 
this molecule keeps the orthorhombic space group (Pbcn, 60) during the spin 
transition, the cell parameters remaining close and the molecular packing the 
same. In the case of [Fe(btz) 2 (NCS)2], the SCO transition is not accompanied 
by a structural transition [11J. Although molecular arrangement in each of 
the two complexes is comparable: layers parallel with the a, b plane (xOy 
basal plane). The intermolecular interactions at ambient temperature and 
their reinforcement at low temperature seem primarily of intra-layer type 
in [Fe(phen)2(NCS)2] and inter-layer in [Fe(btz)2(NCS)2]. This underlines a 
coupling of different nature between close molecules which could be at the 
origin of cooperative effects [H]. Contrary to [Fe(phen) 2 (NCS)2], the SCO 
transition in [Fe(btz) 2 (NCS) 2 ] is gradual and T^ is about 215 K [T2] . 
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2 Theoretical and computational frameworks 

Although the use of Hartree-Fock (HF) approach has been shown to provide 
a good description of the molecular orbital and chemical bonding properties 
mainly in organic chemistry, it becomes well established that calling for the 
density functional theory (DFT) framework [13J brings far more accurate 
results regarding the energetics and related properties. This is because the 
compulsory exchange and correlation (XC) effects are equally treated, albeit 
at a local level, within DFT while only exchange is well accounted for in 
HF although in a better way (exact exchange |14j ) than in DFT. Taking 
the best out of each one of the two approaches led to improvements in ab 
initio molecular calculations with the so called "hybrid functionals" . They 
consist of mixing exact HF exchange, ex. following Becke [H] and DFT 
based correlation, ex. following Lee, Yang and Parr, i.e., the so-called LYP 
correlation [15], with proportions that help to reproduce molecular properties 
of several systems. The calculations have been performed using Gaussian03 
package [To] . 

In the solid state the molecule is built within the framework of an ex- 
tended solid, according to the X-ray diffraction crystal structure determina- 
tion [HIE]- This allows for a better account of possible cooperative effects and 
for a description of the electronic band structure of the whole system. For 
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this purpose we use the all electron augmented spherical wave method (ASW) 
[T7j. Beside its use of DFT, ASW method is based on the atomic sphere ap- 
proximation (ASA), a special form of muffin-tin approximation which consists 
in dividing the cell volume into atomic spheres whose total volume is equal 
to the cell volume. The calculation being carried out in the atomic spheres 
space, empty spheres (pseudo-atoms) need to be introduced -without symme- 
try breaking- in open (low compactness) structures such as that of the system 
studied here. Further, empty spheres allow for the iono-covalent characterist- 
cis of the system to be accounted for by receiving charges from actual atomic 
species. The use of this method in molecular systems was formerly validated 
by Eyert et al. [T5]. In the context of this study, ASW-LDA was used in 
order to obtain a description of the electronic structure with the partial, 
site projected, density of states. Further the chemical interaction is assessed 
especially within the coordination sphere of divalent iron thanks to Ecov 
(covalence energy), which allows to get two-body chemical bond character- 
istics [19]. Negative, positive and zero Ecov point to bonding, anti-bonding 
and non-bonding interactions respectively. 
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3 Results and discussion for the isolated 
molecule 

Starting geometries for HS and LS configuration molecules were extracted 
from the Crystallographic Information File (CIF) files [8]. After testing dif- 
ferent functionals and basis sets through preliminary trial computations, we 
found most accurate results for the energy and vibrational spectra in LS 
and HS states with hybrid functional B3LYP and large and complete enough 
basis sets such as so-called 6 — 31G and 6 — 31 1G, as well as LanL2DZ 
basis sets. The latter include double-C with the Los Alamos effective core 
potential for Fe, S and the Dunning-Huzinaga all-electron double-^ basis set 
with polarization functions for the H, C, and N atoms (( is the exponent 
in the Gaussian type orbitals GTO) [201 EI]- Firstly the geometry of the 
molecule is optimized in both spin states and the results tested against the 
experimental distances within the coordination sphere of iron (cf. fig. [T]). 
Table 1 gives the obtained distances and angles for LS and HS configuration 
from B3LYP / LanL2DZ calculations; we note however that close values were 
equally obtained with 6-311G basis set. Although the final Fe-N distances 
do not translate a perfect octahedral O^-like environment, we find in the 
average d L s(Fe - N) ~ 1.990 A, smaller than d HS (Fe - N) ~ 2.186 A. This 
is an expected result due to the increase of the volume of divalent Fe coor- 
dination sphere arising from the occupation of anti-bonding e* states with 
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two unpaired electrons: LS (t| 9 e°) HS(tf g e^). These average values agree 
with theoretical ones for other similar SCO complexes in the literature [U] 
where BP86 hybrid functional and triple-^ basis set were used. Similarly, 
an agreement between computed values of distances versus experiment was 
also found. However some discrepancy was observed for the magnitudes of 
Fe -IV - C and (SCN) -Fe- (NCS) angles, not only with LanL2DZ 
basis set but with all other basis sets used. From the geometry optimized 
configurations with B3LYP / LanLIDZ , the infra-red (IR) and Raman the- 
oretical frequency spectra were computed for the purpose of providing a 
spectroscopic signature of the system in its two spin states. The differences 
that can be observed for the positions as well as for the relative intensities 
between IR and Raman spectra in both HS and LS states point to the fact 
that there are active modes in IR which are not so in Raman and vice versa 
(cf. figs. Mversus\3§. 

From this we have spectrocopic signatures for the systems not only at the 
level of the spin states (LS/HS) but for the IR/Raman spectroscopies as well. 
For instance this should be particularly helpful when micro-Raman charater- 
izations are carried out. Firstly we discuss the domain of low frequencies 
in wave numbers: cm -1 . At v < 500 cm -1 the Fe-N elongation modes are 
present. In the case of LS, a peak at 532.1 cm" 1 could be detected and as- 
signed to Fe-N(jvcs) elongation which is mainly active in IR. This magnitude 
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comes close to experiment; i.e., 532.6 cm -1 [HI E]). For the HS case this 
same elongation is found at v =240.7 cm -1 (exp. 249.0 cm^^E]). Fe-N( blz ) 
HS/LS elongation values are 221/367 cm -1 ; they follow the same trend as 
above for the magnitudes however we have no experimental data available 
to confront with. The other relevant frequencies characterizing the HS/LS 
varieties which could serve as a signature, are those of NCS~ ligand; their 
magnitudes and attributions are given together with experimental results in 
table 2 with a good agreement. Values for LS state are systematically higher 
than for HS which can be understood if one considers elementary citeria 
such as the average distances which are shorter for the LS complex (cf. table 
1). The frequency of the electromagnetic wave which induces the vibration 
of elongation is given by the relation ^—^y^ where k is the constant of 
bonding strength (considered here as a spring), proportional to the binding 
energy and \x the reduced mass of the two atoms connected by this bond. 
The distances in the case of LS are smaller than in HS case; thus the ratio 
k('fj-) = 0.962, has a magnitude lower than 1. 

At this level of development of the results, it becomes relevant to dis- 
cuss the energetics, A ls /hsi between the two spin states with the functional 
and basis set utilized. The average value sobtained for A HS _ LS is ~ -2.5 
kJ/mol favoring the high spin case. We note that the A LS / HS is governed 
by the amount of exchange admixture in the energy functional [221 [22] • A 
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main result was that adjustment of the exact exchange admixture is neces- 
sary since pure density functionals tend to favor low-spin states, while hybrid 
functionals like B3LYP usually give the high-spin state as the lowest-energy 
state if the low-spin/high-spin energy splitting Als/hs is not too large [21]. 
Our results are in agreement with those of Paulsen, Trautwein and co-works 
[5] who have conducted an extensive study on several spin-crossover com- 
plexes using different basis sets and density functionals. Their results show 
that the B3LYP favors, in general, the high-spin state while all pure density 
functionals favor the low-spin state. Interestingly we shall show that the 
energetics extracted from the solid state calculations (section 4) will exhibit 
this tendency. 

Further we have tried to reproduce the T| characteristic temperature 
from the basis sets that were used. At the transition which is mainly driven 
by vibration energies, equal proportions of HS and LS populations are co- 
existing and one defines from AG = AH -TAS; where G, H and S are 
the well known letters indicating free energy, enthalpy and entropy. At this 
equilibrium AG = 0, and T leads to T~ defined as The electronic 

contribution to the entropy arises from the HS configuration since AS e i = 
R(Ln(2s+l)#s -Ln(2s+l)is; we signal total spins by script s in order to 
differenciate from entropy S; further R = 8.3 J [UJ. Since all spins are paired 
in the t 2g manifold, t^, S= while it s= 2 for HS. Then S to t = Shs — 
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R Ln(2 s^s+l) = 13.4 J. mol -1 ^ -1 . The computed value we obtain is 
13.38 J. mol -1 /^ 1 which comes very close to the formal one. Adding S e ; to 
^vibrational leads to compute T| ~ 100K which is about half the experimen- 
tal one (215K). Such a discrepancy can be well obtained with calculations 
on the single molecule such as [Fe(phen)2(NCS)2] for which a magnitude of 
1530K (T|(exp.)=176 K [7]) was obtained in the literature with BP86 hybrid 
functional |22j . 

One can interpret the reduced T| value from the fact that we have a 
calculated AH = 6.4 kJ/mol while the magnitude we extract from the ex- 
perimental magnetic susceptibility, \m T = f(T) curve [TT], amounts to 
24 kJ/mol. Further AS is smaller than that obtained by the thermody- 
namic relation AG = AH -TAS which gives AS— 112 joules while it is 66 
joules whith the 6-31g basis-set calculations. This raises modeling problems 
in reproducing experimental magnitudes of T| which can have several ori- 
gins, one being the cooperative effects which call for the use of neighbouring 
adjacent molecular species beside the single isolated molecule. We expect 
improvements of T~ provided such neighbouring molecules be adjoined to 
the complex. The calculations are underway. 
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4 Calculations for the extended solid 

To gain insight into the role of each atomic species in the electronic band 
structure and the chemical bonding, we carried out all electrons, spin de- 
generate, calculations for the tetra-molecular [Fe(btz)2(NCS)2] system in the 
orthorhombic Pbcn experimental structural parameters [HI E] . Our calcula- 
tions show large charge transfer, translating ionic behaviour, from the atomic 
species Fe, N, C and S to the empty spheres (ES). Remember that ES are 
pseudo-atoms with zero-core, introduced as stated above, for an account of 
the large voids within the open structure of the molecular system. The plots 
of the projected density of states (PDOS) for the constituent atom sites are 
presented in figs. 4 and 5 for LS and HS configurations respectively. Along 
the the energies are considered with respect to the top of the valence 

band VB (Ey), not the Fermi energy (Ep), because the system is close to an 
insulator, i.e. there is a DOS minimum at this energy boundary, although the 
gap can be hardly seen. This drawback can be attributed to the use of plain 
LDA functional based on the homogeneous electron gas. It is worth saying 
here that the insulating character is likelier to be obtained if LDA+U heavier 
calculations were carried out; U is the on-site Coulomb repulsion parameter 
which enhances the correlation. However the purpose of our work is to un- 
derstand the different behaviours of our system in its two spin states. This 
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is done by discussing the PDOS which are most significant in the neighbour- 
hood of Ey. For this reason we consider a narrow energy window, -6,6 eV 
range. There can be seen a dominant presence of Fe PDOS (red solid lines) 
on both sides of Ey. Its chemical bonding with other constituents species, 
especially N, is observed in the lower part of the VB, see next section. The 
original interesting feature which can be traced out is relevant to the intrinsic 
properties of the complex system in LS and HS states. The major different 
feature between figs. 4 and 5 is relevant to the splitting of Fe(d) states ac- 
cording to the crystal field expected for LS and HS, i.e. the empty e 9 -like 
manifolds PDOS are high in energy within the conduction band (CB) (fig. 

4) in LS with a separation of ~2.5 eV or 242 kJ/mol. which comes close to 
typical crystal field splitting magnitude A in octahedral-like environments. 
On the contrary the partially filled HS e 9 -like manifold is close to Ey (fig. 

5) . 

The variational energy difference between HS and LS ASW-LDA calcula- 
tions carried out with the same Brillouin zone k mesh precision amounts to 
AEls-hs = -8374.6 J for 4 formula units (fu), i.e., -2.09 kJ/fu. While this 
favors the LS configuration energetically, contrary to all molecular calcula- 
tions above, it is interesting to note a close magnitude to the average value 
obtained from single molecule calculations. 

Regarding the chemical bonding within the complex system one mainly 
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needs to examine the Fe-N bond. Considering the LS configuration (the same 
reasoning holds for the HS one), table 1 gives the shorter d,pe-Arcs = l-979 
A distances than d Fe _ N ( btz ) (1.987 and 2.005 A) so one expects relatively 
stronger bonding for the former. The ECOV plots shown in fig. 6 exhibit a 
mainly bonding (negative magnitude) Fe-N character within the VB. How- 
ever a striking feature of a largely dominant Fe-NArcs interaction is mainly 
present. While one can assign a certain role to the smaller d.Fe-Ncs as with 
respect to the two other distances, leading to a stronger Fe-N bong between 
Fe and NCS ligand, one needs to say that the major bonding is due to the 
charged ligands only. This gives the following picture "Fe 2+ (NCS") 2 " (btz) 2 , 
i.e. with the btz ligand playing a secondary role in the coordination sphere 
of Fe 2+ . 

5 Conclusion 

In this work, we presented a molecular calculation for the [Fe(btz) 2 (NCS) 2 ] 
complex. The change of metal-ligand bond lengths upon spin-flip from singlet 
to quintet which is rationalized in terms of a one-particle picture as an occu- 
pation of molecular orbitals with anti-bonding character, is well-reproduced 
in the calculations at the level of the single molecule. The calculated IR and 
Raman absorption spectra for the low-spin and high-spin, we obtained with 
the LanL2DZ frequency of antisymmetric and symmetric NCS- were very 
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close to those obtained by experiments. We also carried out original calcula- 
tions for the solid states accounting for the whole crystalline tetra-molecular 
complex. With the all-electrons ASW-DFT method we obtained the PDOS 
in the two states and the chemical bonding characteristics. The first show in 
the close neighbourhood of the top of the valence band the change from LS 
to HS behaviour within the Fe coordination sphere for the crystal field. The 
chemical bonding shows the original feature of the predominant Fe 2+ -NCS~ 
bond with respect to much lower magnitude Fe-N(btz) crown. 
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LS 


LS exp.[8, 9J 


HS 


HS cxp.[8,E] 


d-Fe-NCS 


1.979 


1.947 


2.085 


2.066 


&Fe-N{btz) 


1.987 


1.967 


2.267 


2.177 


&Fe-N{btz) 


2.005 


1.972 


2.207 


2.166 


a(Fe-N-C)° 


155.20 


169.1 


145.30 


167.80 


a(NCS) -Fe- (NCS)° 


96.132 


86.06 


107.72 


91.53 


a{N htz ) -Fe- {N btz )° 


80.0 


80.0 


73.92 


74.7 


(dpe-N) 


1.970 


1.979 


2.186 


2.136 



Table 1: Distances in A and angles compared to experiment obtained by 
using LanL2DZ basis set and B3LYP hybrid functional for the molecule 
geometry optimisation in LS and HS [Fe(btz) 2 (NCS)2] complex. Average 
(dFe-N) distances are provided for sake of discussion (see text). 



LS IR LS Raman 


HS IR HS Raman 


2103.89 (2100) 2103.89 (2107) 
2113.51 (2110) 2113.51 (2110) 


2070.71 (2060) 2070.71 (2062) 
2077.78 (2070) 2077.78 (2072) 



Table 2: Comparison of antisymmetric and symmetric N-C-S frequencies in 
cm -1 between LS and HS [Fe(btz) 2 (NCS) 2 ] complex. Experimental values 
from [31 E] are given between brackets. 
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Figure 1: (Color online). Sketch of the molecular structure of SCO complex 
[Fe(btz) 2 (NCS) 2 ] with characteritic distances (A) and angles shown here for 
the low spin LS case. 
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Figure 2: IR spectra of HS (in the top) and LS (in the down) of 
[Fe(btz)2(NCS)2] complex with the LanL2DZ basis-set. 




Figure 3: Raman spectra of HS (in the top) and LS (in the down) of 
[Fe(btz) 2 (NCS) 2 ] complex with the LanL2DZ basis-set. 
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Figure 4: (Color online). Site projected density of states (PDOS) of LS 
[Fe(btz) 2 (NCS)2]. Energy reference is with respect to the top of the valence 
band (Ey). 
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Figure 5: (Color online). Site projected density of states (PDOS) of 
[Fe(btz)2(NCS)2] in the high spin (HS) case. Energy reference is with re- 
spect to the top of the valence band (Ey). 
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Figure 6: (Color online). Chemical bonding following ECOV criterion (arbi- 
trary units) for the different Fe-N interations within [Fe(btz)2(NCS)2] system 
in the low spin (LS) case. Energy reference is with respect to the top of the 
valence band (Ey). 



